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SUBJECT} A PROGRAM TO SOLVE THE CLASSICAL TRANSPORTATION PROBLEM 

A program has been written to solve the classical transporta* 

1 
tion problem on the Whirlwind computer using the stepping stone method,, 

The transportation problem may be stated as follows s A eom- 

th 
paay operates m plants producing a commodity, the i of which can sup* 

ply S. units of the commodity* The company sells its production t© n 

customers, the j of which desiree D. uaits of the commodity, The 

cost of manufacturing and transporting a unit of the commodity from 



plant i to customer j is C. .« It is desired to find the number of 
units X. . that should be shipped from each plant to each customer, 
order to have the total cost of the operation be a minimum* 



Thus the problem may be stated mathematically ass 




minimize G «* r"".C 4 .X, . 


(1) 


subject to the constraints ^ 




J X iJ- D J 


(2) 


i X ij * S i 


(3) 


X iJ 2 ° 


W 



This is a special ease of the general linear programming problem, 

A set of X.. which satisfies equations 2, 3, and 4 is called 
a feasible solution to the problem* If the set also minimi zes C, it 
is an optimum feasible solution containing, in general, exactly m+n-1 



non-zero X. .- 



The proeess of solution consists of s 

1) generating a feasible solution having exactly m+n-1 
non-zero X. . • 
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2) finding a zero X.. which, if allowed to be positive, 
would yield a dedrease in G 6 In the process of in- 
creasing this X.., one of the non-zero X.. must feo 
to zero in order that equations 2 and 3 remain satis- 
fied* Thus the new feasible solution will again have 
exactly m+n-1 non-zero X. .» 
Step 2 is repeated until there are no zero X.. which can be' 
changed so as to reduce C« The set of X.. is then an optimum feasible 
solution. 

The .Whirlwind program is unique in that a logical search pro- 
cedure is used in entering the new element into the feasible solution 
whereas other transportation codes use auxiliary tables through which 
many passes must be made for each step towards the solution., Also, in 
finding which zero X.. should be made positive, the usual procedure is 
to search thrbugh all of the cost data to find the X. . which could b© 
mad© non-zero to the greatest advantage* The Whirlwind program takes 
the first X,. . found that would lower the total cost if made noa-zer© 
This avoids much of the slow and time consuming process of searching 
through cost data at the expense of a greater number of iterations* 

Because of these features this program is considerably faster 

2„ 
than other procedures A number of problems with 9 plants and 69 

customers have been solved in about 1 5 minuses of computer time each. 
The average number of iterations was 150 « Also, seven problems with 
m ■ 60, n * 291, have been solved requiring about 45 minutes and 1800 
iterations each. 

The program will handle problems with m £ 127 , m+n £ 401* 
In many practical problems not all mn possible G . u are significant* 
Ho error will result if the costs for shipping routes that are known 
to be absurd or impractical are assumed to be infinite* Th© program 
has therefore been arranged to consider as infinite all costs not spe- 
cified in the data. A considerable saving in storage is ma4© this way* 
Problems with up to 10,000 silgnificant costs can be solved by th© pro- 
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The appendix to this report contains detailed instructions 
for the preparation of data tapes for the transportation routine and 
directions for its use„ 



References g 

1* A. Gharnes and W«W. Cooper, "The Stepping Stone Method of Uxplalniag 
Linear Programming Calculations in Transportation Problems," 
Management Science , vol* 1, no« l t October 1954* 

2. For example, see C W Swift and S* Poley, "The Transportation 
Problem," IBM Technical Newsletter No. 10, October 1955. 



DGL-140 Page 4 

Appendix 

Content a s 

A« Data Tape Preparation Page 4 

Bo Output Format Page 8 

C» Performance Requests Page 9 

A« Data Tape Preparation 

The data that must be supplied to the routine to solve a 

problem includes? 

a« the size of the problem (m and n) 

b the cost data G. . 

ij 
Co Supply capacities S. 

d» Customer demands D. 

The data tape for a sample problem is given in Fig* A-l 
Specifications for the data tape are given below* The symbol -^ denotes 
a tabulation, -) denotes a carriage return a The important parts of the 
data tape are discussed in their order of appearances 

lo 

"f c TAPE xpxxx-xxtqcx" 

Tape number and title in conventional GS II format, 
2o 

"32 a .booed" Output Code 

a ■ Cost data C.. are supplied on this 

tape 
a » 1 Cost data of the previous problem 

are to be usedo 
b = Supplies and Demands are given on 
this tape» The computer will stop 
after reading the data tape« 
b ■ 4 Supplies and Demands are not given 
on this tape c The computer will, 
stop after reading the data tape* 
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c ■ Output will be punched on paper 

tape* 
e ■ 1 Output will be print ed* 
d ■ Final basis only* 
d ■ 2 Initial and final bases* 
d ■ § Initial and final bases and step 

data* 
d ■ 4 Final basis and alternate solution 

element So 
d ■ 6 Initial and final bases and alter- 
nate solution elements, 
d ■ 7 Initial and final bases, step data 
and alternate solution elements* 

3. 

"STABT AT 56" 

(One and only one carriage return!) 



4. 



5. 



"Sample Problem Ho* 000" 

Title of problem, (This title will appear 

on the results*) The last character of the title line 

must be a carriage return* Any flexowriter characters 

except a carriage retum»may be used elsewhere in the 

title* 
* 

*+ a + n * ■ 
Size of Problem 



^Between the asterisks the following characters are ignored 
by the computers », + , eolor shift, upper case , lower case, space* 
back space, stop character and delete. All others are illegal except 
for the digits © and 9 and ci|aracters specifically mentioned* Also 
a word enclosed in parenthesis will be ignored* i.e* "(costs for plant 
2)*" Any flexowriter characters may be used in the word except tabula- 
tion or carriage return* The "comment" word must be followed by a 
tabulation or carriage return,, 
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m » Number of supplies ^.127 
a ■ Number of demands ^z. 401-m 
6. 

Cost data C. . if included im this tape. 
Ther&'^e f ourolegal word t#pes which' must be . separated 
by -at least one tabulation- or carriage return. 

a) M + .08964" 

Cost data word« Plus sign and 
decimal point are ignored. The decimal 
point is assumed to be at the extreme 
left „ hence ^numbers must be less than 
+1.0. The numbers 9 which are assumed 
to be positive, must have no more than 
9 digits including zeros* The first 
cost data word is assumed to be C__ un- 
less a j assignment word proceeds* 

b) "j m 28" 

j assignment word* The next 
cost data word will be assumed to be 
for column 28 in this row. The j value 
specified must be less than or equal t© 
n and must not be zero* The equals 
sign may be omitted, 
e) »| w 

Row termination word. Denotes 
the end of this row. The next cost 
data word will be assumed to be for 
column one in the next row unless a j 
assignment word intervenes « 
d) "x" 

Termination word. Denotes the 
end of the last row of cost data. 
(The vertical bar may not be used t© 
terminate the }.ast row.) 
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For any row of C. , twice the number of 
cost <|ata words plus the number of j assignment words 
must be less than or equal to 2n. The total amount of 
eost data must be less than 13,888 registers, that is? 

2c + m + j £t 13,888 [ (A-l) 

Here, c is the number of cost data words, j is the number 
of j assignment words and m is the number of rows in the 
0, . matrix* Costs which are omitted from the data tape 
are assumed to be "infinite" ( +1*0 ) by the program© 
7. 

Supplies and Demands, S ± and Dj if in-. 
eludeip. in this tape* There are three legal word types 
which must be separated by at lea&b one tabulation or car- 
riage return* 

a) Supply or Demand wordo -Same specifi- 
cations as for cost data word above* 
The first m- such words will be the sup- 
plies, S.*«oS « The next n such words 

l m 

will be the demands D, *„*D - 

i n 

b) n |" 

Denotes the end of the supply data 
words* 

c) "x" 

Denotes the end of the demand data 
words « 

* (See note, Page -5) 
8* 
"fen 
STAp AT j w 

End of tape* (The *fe* may be the tiegin- 
ning of another data tape which will be processed automa~ 
tienlly after the first tape*) There must be n© flexo- 
wr^ter character other than feed out between the n x n 
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terminating the data and the "fc"* About one inch of 
feed out is necessary ,. One carriage return must follow 
"3TART AT 0"* 

B* Output Format 

Fig* A-2 shows the output for the sample problem* The 
"Total supply" and "Total flemand" are given as a check on the data. 
They should not be expected to agree exactly because of error in con* 
verting from decimal to binary* 

The "Initial basis table" gives the location and value of 

the X. . in the initial feasible solution. The U's and T»s are extra 

variables that are useful in obtaining the solution* Their values 

are such that TL*7. * 0. . for all combinations of i and j for which 

X ' is n< 

ij 
by rows. 



X 'is non-zero in a feasible solution,, They are arranged in order 



The step data gives for each step the row and column of the 
noa-zer© X. . dropped and the non-zer© X. . added to the feasible s©lifo* 
tioa, the value of the new X. . and the incremental negative cost of 
making the change. 

^h© "Initial cost" is the sum of \JL** *©* "^he initial > 
feasible sdiutioa. The "Total decrease" is the sum of the decreases 
in cost at ^ach step* The "Final cost" is obtained by subtraction© 
The "Summed cost" is given for comparison and is the sum of GL .X.. for 
the optimum solution. This will not be exactly equa^. to the final 
cost because of round-off error in multiplication. 

1*ae number of "Iterations" is the number of steps required 
t© ©btain the optimum solution from the initial solution,, The number 
of "Passes" is the number of times the whole set of cost data had t© 
be searched to find all of the zero X.. that could b$ brought int© the 
solution profitably* 

The "Alternate elements" are those shipping routes which, if 
used, would not increase the total cost* 
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C . Pe rf ormance Re que s t s : Operator Instructions 

Once data tapes have been prepared, the problem solutions 
may be obtained by following the instructions below which are given in 
the standard form for the Whirlwind machine; 

1) si 1 switch off 

2) e, fb 219 - 68 - 2000, ri Program tape 

3) fc xxx-xx-xxxx, ri *) 

----- y Data tapes 

fc xxx-xx-xx, ri J 

As many data tapes may be run as desired on a single read-in 
of the program tape. One read-in of a set of cost data may be used 
for any number of solutions. However, supply and demand data are des- 
troyed by the routine and, hence, must be given for each run. 

If the quantity 2c + m * jvin Equation A-l of section A on 
data tape preparation is greater than 13,888 but less than 20,480 the 
problem can still be handled by transferring the program to the buffer 
drum. The operating procedure is as follows: 

1) si 1 switch off 

2) e, fb 219 - 68 - 2000, ri Program tape 

3) fb 219 - 68 - 3000, ri ^ Tape to transfer 

program to buffer 
drum, 

4) fc xxx-xx-x C Data tapes 

If this procedure is used, parts 1, 2 9 3 and 8 of the data 
tape must be in 556 binary form. Otherwise the transportation routine 
on the buffer drum will be destroyed by the CS II conversion program. 
Performance Requests : Input-Output Equipment 

The amount of input-output equipment used by the program is 
as. follows: ■ , . ' ■ 

Direct typewriter - two lines for each set of supplies and 

demands read in. 

Auxiliary Storage Drum - groups 0-10. 
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Magnetic Tape - unit #3. The amount of magnetic tape used 
may be found from the following equation. 
Number of feet ■ N 

^2l L* 1 + a)U + a)(!8 + 30C) + b(25 + 2(£)@ + c(9y) + 50o] 

(A-2) 
Where s a * 1 if initial basis is requested 
b * 1 if step data is requested 
c * 1 if alternate elements are requested 
OC *= Maximum number of digits to the right of deci- 
mal point in the data 
{3 « Number of iterations. 
"V"« Number of alternate elements* 

Performance Requests ; Time Requirements 

The approximate time required to solve a given problem may 
be estimated from Fig. A-3. About one minute is required to read in 
the program tape. 

Performance Requests ; Alarms 

There are three possible printouts on the direct typewriter 
that can result froiai a faulty data tape. Their possible causes are 
listed below : s 

1. "wrong count" 

a. To© much cost data was given in this 
row of the G. . matrix. 

b. Exactly m rows were not given for the 
cost matrix. 

c. Exactly m supplies and exactly n de«* 
mauds were not given. 

2. "illegal character" 

a. The output code specified in part 2 of 
the data tape was illegal* 

b* As. illegal flexowriter character was 
used improperly* 
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3. "number too long" 

a» More than 3 digits appeared in ra or a. 

b, A j assignment value greater tham m 
was used* 

c« More than 9 digits appeared in a cost 
or a supply or demand word* 
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FigUTe A-l 



fc TAPS 219 - 68 - 000 DENNIS 
3^ 0.00007 



START AT 56 

Whirlwind Transportation Routine - Sample Problem 
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Whirlwind Transportation Routine 



Total supply 
Total demand 



.7599 
.7599 



Initial basis table 



.0399 
.0000 
.1899 
.0599 



Sample Problem 



Fig. A-2 



5 


2 


.0799 


1 


5 


.0599 


4 


8 


.0399 


4 


4 


.1399 



X 

.0499 
.0099 
.0699 
.0199 



.0000 



.9915 



.9875 



.9955 



.9985 



.9885 
- .9875 




- .9915 
+ .0080 


Step data 
11 Jl 

3 4 
2 3 

1 5 

2 7 
5 6 
2 1 
1 3 

4 4 

4 5 
1 2 
1 3 

5 1 


12 J2 

3 5 

4 3 

f 5 

4 7 

5 1 
1 2 
1 3 
1 4 

1 5 

2 1 
2 3 
5 6 


Initial cost 
Total decrease 
Final cost 


.076005815 
.070595821 
.005409994 


Final basis 

1 J 
2 
2 
1 
3 


table 

1 
4 
7 
5 


.0399 
.0700 

.1199 
.0599 


Summed cost 




.005409994 


U 

+ .97^5 




+ .9755 


V 

- .9725 

- .9625 




.9665 
.9665 


Iterations 
Passes 


3 


12 


Alternate elements 

1 J ■ 
3 * 


1 



.9835 
.9905 



X 

.0599 
.0499 

.0000 

.0099 
.0199 
.0199 
.0499 

.0400 
.0000 

.0199 

.0499 
.0199 



.9715 



.9675 
.9655 



.9805 



.9905 
.0030 
.9965 
.9995 
.0049 
.0040 
.0020 
.0090 
.0030 
.0049 
.0039 
.0040 



.OO69 



x 

.0799 
.0000 

.0399 
.1299 



.9705. 



.9645 



2 


i 
3 


X 

.0499 


4 


6 


.0099 


4 


7 


.1399 


5 


6 


.0199 



.9735 



.9675 
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Time Required for Solution 
Fig. A-3 
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